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Abstract 

The nuclear spin-dependent parity nonconserving (PNC) interaction arising from a combination 
of the hyperfine interaction and the coherent, spin- independent, PNC interaction from Z exchange 
is evaluated using many-body perturbation theory. For the 6,Si/2 — 7si/2 transition in ^'^^Cs, we 
obtain a result that is about 40% smaller than that found previously by Bouchiat and Piketty 
[Phys. Lett. B 269, 195 (1991)]. Applying this result to ^^SCs, leads to an increase in the experi- 
mental value of nuclear anapole moment and exacerbates differences between constraints on PNC 
meson coupling constants obtained from the Cs anapole moment and those obtained from other 
nuclear parity violating experiments. Nuclear spin-dependent PNC dipole matrix elements, includ- 
ing contributions from the combined weak-hyperfine interaction, are also given for the 7si/2~ 85^/2 
transition in ^^^Fr and for transitions between ground-state hyperfine levels in K, Rb, Cs, Ba^, 
Au, Tl, Fr, and Ra+. 
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I. INTRODUCTION 



The precise measurements of the 6s[F = 4] — 7s[F = 3] and 6s[-F = 3l — 7s[-F = 4] parity 
nonconserving (PNC) dipole matrix elements in ^^'^Cs by Wood et al. [l| lead to a value of 
the weak charge^ Q'^^{^^^Cs) = —72.73(46) that is in agreement with the Standard Model 
value QiY^C^^Cs) = —73.09(3) 11]. These measurements also lead to an experimental value 
of the much smaller contribution from the nuclear spin-dependent PNC interaction that is 
accurate to about 15%. This spin- dependent contribution has three distinct sources: the nu- 



clear anapole moment 



mm, the Z exchange interaction from nucleon axial-vector (AnVe) 



currents, and the combined action of the hyperfine interaction and the spin-independent 



Z exchange interaction from nucleon vector (V^Ag) currents IJ, |l5|. Of these three, the 
anapole contribution dominates. The contributions from the anapole and nuclear axial- 
vector current are 

if» = ^«:,a.Jp(r), (1) 

where G is the universal weak coupling constant, / is the nuclear spin, and p(r) is a nor- 
malized nuclear density function. The subscript i of the dimensionless constants Hi takes 
the values i = a for the anapole contribution and i = 2 for the axial- vector contribution. In 



Refs. 



15l |. the hyperfine-vector current contribution was also reduced to the form given 
in Eq. (Q) with a corresponding dimensionless constant k^- 

To extract the anapole contribution Ha from experiment, it is necessary to know the 
corresponding spin-dependent PNC amplitude calculated with Kj(i — 1, clS well as the two 
contributions from the axial-vector and weak-hyperfine interference terms quantified by k,2 
and K]^ Th e spin-dependent PNC amplitude was calculated in various approximations in 



ihf. ±ue s 



Refs. I3I lid, Il7|. Nuclear shell-model values of ko for ^^^Cs and ^os^jj ^ere obtained in 

by Flambaum and Khriplovich [1^ and values of were later determined for various cases 
of experimental interest by Bouchiat and Piketty (iq . 



Recently, Haxton and Wieman used the values of K2 and hm determined as described 



^ The experimental value includes a net correction of -1.1% to the theoretical PNC amplitude y,y,y| from 
the Breit interaction g aZ vertex corrections flH, Coulomb-field vacuum polarization % and nuclear 
skin eff ects 

m. 

^ In Ref . [l9j , K2 and Khf are designated by kzo £^nd kq„ , respectively. 
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above to extract values of Ka from PNC measurements in ^^^Cs fl^. The resultin g an apole 
moments were, in turn, used to place constraints on PNC meson coupling constants 20| . The 
constraints obtained from the Cs experiment were found to be inconsistent with constraints 
from other nuclear PNC measurements, which favor a smaller value of the ^'^'^Cs anapole 
moment. 

Motivated by this disagreement, we are led to re-examine the combined hyperfine-weak 
interaction. We find that the contribution of this term to the PNC dipole matrix element at 
the Dirac-Hartree-Fock (DHF) level can be approximated by a spin-dependent interaction 
of the type given in Eq. (0); however, such an approximation is not justified in correlated 
calculations, since contributions from Eq. are very sensitive to correlations, whereas 
contributions from the combined hyperfine-weak interaction are relatively insensitive to cor- 
relation corrections. We do find, nevertheless, that even in correlated calculations there is a 
rough proportionality between contributions from the combined interaction and those from 
the interaction given in Eq. that is independent of hyperfine state, and we use this fact 
to define "effective" values of the coupling strength /thf for cases of potential experimental 
interest. For the 6s — 7s transition in Cs, our effective value of km is about 40% smaller than 



the value from 15|. Interestingly, for this case our final correlated value of is quite close 
to the value predicted by the formula derived in Other things being unchanged, this 
decrease in the size of K^i leads to an increase in the size of and, correspondingly, in the Cs 
anapole moment; consequently, increasing the inconsistencies between various experimental 
constraints on PNC meson couphng constants described in jioj ]. 



II. METHOD 



We write the spin-dependent PNC correction to the reduced electric-dipole matrix ele- 
ment {wFf\\z\\vFi) as the sum of three terms: 

{wFF\\z\\vFj)f^c = >^a {wFf\\z\\vFi)^'''^ + K2 {w Ff\\z\\vFi)'^^^ + {wFf\\z\\vFiY'''\ (2) 

where indices (a), (2), and (hf) correspond to the anapole, axial- vector, and weak- hyperfine 
interference, respectively. Since the anapole and axial-vector contributions both take 
the form given in Eq. (^, we can introduce {wFf\\z\\vF]:Y'^'°'^ = {wFf\\z\\vFi)^'^'> = 
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{wFp\\z\\vFj)^'^\ We then define Khf as the ratio 

_ {wFp\\z\\vFj)('^'^ 
"""^ {wFF\\z\\vFr)(^^-)- 

We expect, and indeed find, that km depends on the initial and final hyperfine levels. For 
cases considered here, however, the dependence of on the hyperfine levels Fj and Fp 
is weak and we may treat as constant to some level of accuracy. We write the expres- 
sion for the total spin- dependent PNC contribution to the electric-dipole matrix element 

{wFf\\z\\vFi) as 

{wFp\\z\\vFjYp%c = i^a + K2 + Km){wFf\\z\\vFi)^'^''\ (4) 

and define 

K = K,a + K,2 + Khf. (5) 

In this work, we calculate both {wFp\\z\\vFj)^^^ and {wFp\\z\\vFi)^'^'°'^ and, therefore, de- 
termine the state-dependent values of Khf. 

The hyperfine interaction Hamiltonian is written 

^ ' A-K J \r -xY' 

where M{x) is the magnetization density, which is related to the nuclear moment /x/ by 

Here = \e\h/2Mp is the nuclear magneton. 

The dominant, spin- independent, part of the weak interaction is 

where Qw is the conserved weak charge of the nucleus, given at tree level in terms of the 
neutron number A^, the proton number Z, and the Weinberg angle 6^ by Qw = Z(l — 
4 sin^ 6\y) — N, and p(r) is a nucleon distribution function. In our numerical calculations 
of the interference term, we use radiatively corrected values of Qw inferred from ^|. The 
nucleon distribution p(r) is assumed to have the form 

'^'^ = l + expKr-C)/a] 
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where a = 0.523 fm (corresponding to 90%-10% fall-off thickness t = 2.3 fm) and where C 
is inferred from the nuclear charge radii listed in [21]. In the exceptional case of ^^^Fr, we 
choose C = 6.733 fm, corresponding to the value Rrms = 5.566 fm given in Ref. We 
assume that the radial dependence of the magnetization distribution is identical to that of 
the nucleon distribution. 

As shown in Appendix ^ the dipole matrix element corresponding to the weak-hyperfine 
interference is given by the third-order perturbation theory expression 

{wIFfMf\z\vIFiMi)'^^^^ = 

{w\H^'^^\n){n\H^'^^^m){m\z\v) ^ (wlH^'^^^ln) {n\H'-^^\m){m\z\v) 
"~ {Eu, — Em){Ew — En) ^ {E-uj — Em){Ew — Eri) 

^ {w\H^^'^\m){m\z\n){n\H^'^^^\v) ^ >^ {w\H^^^^\m){m\z\n){n\H^'^^\v) 

J {Ew — Em){Ey — En) ^ [E^ — Em){E^ — En) 

^ {w\z\n){n\H^^^\m){m\H^^^^\v) ^ {w\z\n){n\H^'^^^\m){m\H^'^^\v) 

^ {Ey — Eni){Ey — En) ^ {Ey — Em){Ey — En) 

where we use designations \w) and \v) on the right-hand side for coupled hyperfine states 
\wIFpMp) and \vIFjMj), respectively, and where we designate the energy of state i by Ei. 
Note that the other matrix element {wIFpMp\z\vIFiMj)^'^'"-'> in Eq. is obtained from a 
considerably simpler second-order perturbation theory calculation. 

In Ref. [l5[, terms on the second and fourth lines of Eq. Q were ignored and partial 
sums on the first and third lines, such as 

^ En-Ey ' 

were carried out using free-particle Green's functions and reduced to an effective interaction 

H^^""^ = ^ KM « ■ J p(r) 

of the form given in Eq. (PJ). A similar reduction was made in Ref. Here, we evaluate all 
of the terms in the Eq. Q numerically. One important advantage of this direct numerical 
evaluation is that correlation corrections to {wIFp\\z\\vIFj)^^^ can be determined using 



standard many-body methods. Indeed, we find that Eq. Q is insensitive to correlations at 
the random-phase approximation (RPA) level for most of the cases considered here, whereas 
calculations based on the contracted approximation above are very sensitive to correlation 
corrections! A reduction of Eqs. and Q to reduced matrix elements suitable for numerical 
evaluation is given in Appendix ^ 

III. NUMERICAL RESULTS 

We evaluate the reduced dipole matrix element {wFp\\z\\vFj)^^^' given in (inTill and the 
reduced matrix element {wFf\\z\\vFj)^'^'"''^ given in ()B15|) . and we find that their ratio is 
approximately independent of the angular momentum quantum numbers Fj and Fp for 
transitions between hyperfine levels. 

Let us consider the 6s — 7s transition in Cs. We first evaluate the reduced matrix elements 
in ()B15j) and ()B14|1 at the DHF level of approximation. We solve the DHF equations in a 
finite B-spline basis using the methods described in and use the resulting basis functions 
to evaluate matrix elements and carry out sums over intermediate states. For the case Cs, 
our basis set consists of 100 splines of order 15 for each angular momentum state. The basis 
orbitals are constrained to a cavity of radius 45 a.u.; the cavity radius is modified in other 
atoms to accommodate the initial and final valence orbitals. As a check, we carried out the 
Cs calculations using a cavity of radius 75 a.u. to verify that the results are stable against 
changes in the cavity radius. Results of our DHF calculations for the transitions between 
the possible hyperfine levels are presented in the upper four rows of Table IH We find that 
the ratio km of the {wFf\\z\\vFj)^^^ to {wFf\\z\\vFj)^'^'"'^ matrix element changes from level 
to level by only 2% in the DHF approximation. 

The DHF treatment of PNC in cesium is known to be a rather poor approximation, 
giving a value for the dominant part of the PNC dipole matrix element that is 20% smaller 
than the final correlated value. To obtain a reliable value for the PNC matrix element, one 
must go beyond the DHF approximation and treat correlation corrections. The dominant 
correlation corrections, those associated with core shielding, are obtained in the random- 
phase approximation. Including RPA corrections to both weak-interaction and dipole matrix 
elements gives a value for the dominant PNC dipole matrix element in Cs that is within 2% 
of the final correlated value. 
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TABLE I: We list values of Khf for transitions between hyperfine levels 6s [-F/] — 7s [Fp] in Cs de- 
termined in DHF and RPA approximations. The atomic number is A = 133, the nuclear spin is 
/ = 7/2, the nuclear magnetic moment = 2.5826, the weak charge (including radiative cor- 
rections) is Qw = —73.09(3), and the 50% fall off radius is C = 5.675 fm for both the nuclear 
p{r) and magnetization M(r) distributions; the 10%-90% fall off distance is 2.3 fm. The PNC 
reduced dipole matrix elements {wFf\\z\\vFi)^'^'"'^ are given together with the weak-hyperfine in- 
terference correction to dipole matrix elements {wFp\\z\\vFj)^^^^ ; their ratio is K^f. Numbers in 
square brackets represent powers of 10. 



Type Ff-Fj (7sFp| |z| |6sF/)(2.«) (7sFf| |z| |6sF/)(hf) 



DHF 


3 - 3 


1.908[-12] 


1.193[-14] 


6.251 [-03] 


DHF 


3 - 4 


5.481[-12] 


3.480[-14] 


6. 349 [-03] 


DHF 


4 - 3 


4. 746 [-12] 


3.020[-14] 


6. 364 [-03] 


DHF 


4 - 4 


2.173[-12] 


1.358[-14] 


6.251 [-03] 


RPA 


3 - 3 


2.249[-12] 


1.141[-14] 


5. 076 [-03] 


RPA 


3 - 4 


7.299[-12] 


3.579[-14] 


4.903[-03] 


RPA 


4 - 3 


6.432 [-12] 


3.139[-14] 


4.880[-03] 


RPA 


4 - 4 


2.560[-12] 


1.300[-14] 


5. 076 [-03] 



The RPA matrix elements are calculated as described in Ref. j2^, with the value of 
uj in the RPA equations set to zero. In the last four rows of Table HI we give values 
of {7sFp\\z\\QsFj)^'^'"'^ that include RPA corrections to both dipole and weak-interaction 
operators, and values of {7sFp\\z\\6sFj)^^^'> that include RPA corrections to the dipole, 
weak-interaction, and hyperfine operators. While the RPA values of {TsFfIIzUGsFj)^"^'^^ are 
15-25% larger than the DHF values, the RPA and DHF values of {7sFf\\z\\6sFi)^'^^^ differ 
by only 3-5%. Thus, by contrast to PNC dipole matrix elements induced by the domi- 
nant spin- independent interaction and by the spin-dependent interactions given in Eq. (PJ), 
which are very sensitive to correlation corrections, the third-order matrix elements for the 
combined interaction are relatively insensitive to correlations for the 6s — 7s transition in 
Cs. It should be emphasized that the contraction of operators introduced in j3| is a useful 
approximation at the independent-particle DHF level of approximation; however, when cor- 
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TABLE II: We list values of Khf for transitions between hyperfine levels 8s [Fp] — '''■s[-F/] in Fr 
determined in RPA type calculations. Here, the atomic number A = 211, the nuclear spin is 
/ = 9/2, the nuclear magnetic moment /i/ = 4.00, the weak charge (including radiative correc- 
tions) is Qw = —116.23, and the 50% fall off radius is C = 6.7325 fm for both the nuclear p(r) and 
magnetization M(r) distributions; the 10%-90% fall off distance is 2.3 fm. The PNC reduced dipole 
matrix elements {wFf\\z\\vFj)^'^'°''' are given together with the weak- hyperfine interference correc- 
tion to dipole matrix elements {wFf\\z\\vFj)^^^^ ; their ratio is Khf- Numbers in square brackets 
represent powers of 10. 



Fp-Fi (8sFp||z||7sF/)(2-«) {8sFf\\z\\7sFi)'-^^^ km 

4 - 4 3.092[-ll] 3.472[-13] 1.123[-02] 

4 - 5 1.016[-10] 1.069[-12] 1.053[-02] 

5 - 4 9.224[-ll] 9.645[-13] 1.046[-02] 
5 - 5 3.426[-ll] 3.846[-13] 1.123[-02] 



relation corrections are included, although an approximate proportionality still obtains, the 
proportionality constant depends on correlations; this is a reflection of the fact that there is 
no effective Hamiltonian of form for the combined interaction. 

We include negative-energy contributions when evaluating sums over intermediate 
states in Eqs. ()B14|) and ()B15|) and when calculating RPA matrix elements. We find almost 
no negative-energy correction to {7sFf\\z\\6sFj)^'^'"'\ However, the negative-energy correc- 
tions to (7sFp||2||6sF/)(^^) were found to be large, 22-23% at both the DHF and RPA levels 
of approximation. Since negative-energy contributions are important for accurate calcula- 
tion of {7sFf\\z\\6sFi)^^^\ they are, therefore, important in the evaluation of K^f. Omission 
of negative-energy contributions leads to values of km = 0.0049 in the DHF approximation 
and Km = 0.0038 in the RPA approximation which are about 20% smaller than our final 
values listed in Table IH We note that our final correlated value fortuitously coincides with 
the DHF value without negative energy contributions; the correlation correction decreases 
the value of /thf and the negative energy contribution increases km by approximately the 
same amount. We stress again that these two effects contribute, in fact, to different quan- 
tities, negative energies contribute only to {7sFf\\z\\QsFj)^^^^ and correlation primarily to 
{7sFf\\z\\6sFjY^'^\ 
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We also found that sums in the interference matrix element given in Eq. ()B14|) must 
include the entire set of basis orbitals, in contrast to the sums in Eq. ()B15|) . where omitting 
high-energy orbitals from the basis has very little effect. In other words, the completeness 
of the basis is very important for calculation of the sums in ()B14|) . 

Results of our RPA calculations for the 7s — 8s transitions in Fr are presented in Table ITTl 
The state dependence of km increases to 6-7% in Fr in comparison to Cs, where differences 
in Khf for different transitions were 3-4% in the RPA approximation. As in the case of Cs, 
the largest differences occur between transitions with Fj = Fp and those with Fp ^ Fj; 
there is only 0.7% difference in Khf between the 4-5 and 5-4 transitions. 

n 

For the 4-3 and 3-4 hyperfine transitions in Cs measured by Wood et al. [ij, an effective 
value = 0.0049 can be extracted from the RPA values listed in Table HJ This value is 
about 40% smaller than the value Khf = 0.0078 from but agrees exactly with the value 
obtained earlier by Flambaum and Khriplovich We use our value of to extract a 
value of from the Cs PNC experiment of Wood et al. 



A 



lm{Ei 



PNC) 



= -0.077 ± 0.011 mV/cm, (10) 

J 34-^43 



where (3 is the vector polarizability of the transition, which has been measured in Ref. 
with high accuracy P = 27.02(8) a^. The subscripts 34 and 43 in (fTUj) correspond to FpFj. 

The spin-independent PNC amplitude -Ep^c alkali-metal atoms {jp = ji = 1/2) is 
customarily defined as 

4nc = (JF l/2\z\jj 1/2), (11) 

leading to the following relation between spin-dependent PNC amplitude and the corre- 
sponding spin-dependent reduced matrix elements: 

E^;^^^ = ^{wFp\\z\\vFj)(^-\ (12) 

where the k is defined by Eq. (JSj) and A is an angular coefficient 

A = (^-ly^+Pr+i+i ^6 [Fj] [Fp] ^ \ , (13) 




where [F] = 2F + 1. For the two transitions considered here, A43 = —^34, so we may write 
A 



Im(EpNc) 



34-43 



A43P 



(7sF^||z||6sF,)£"^ + (7sF^||2||6sF,)if^ — (14) 
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TABLE III: Comparison of contributions to spin-dependent PNC in ^^"^Cs obtained by different 
groups. All results are presented in terms of the coefficients k^, K2, km, and their sum k, used in 
the present paper. 



Group 



1^2 



Kn 



Present 

Haxton et al. [18, ^9, J7j 
Flambaum and Murray [28] 
Bouchiat and Piketty [15, 29j 



0.117(16) 
0.112(16)'' 
0.112(16)'^ 



O.OUO'' 
0.0140 
O.Olll'' 
0.0084 



0.0049 
0.0078"^ 
0.0071^ 
0.0078 



0.098(16) 
0.090(16) 
0.092(16)5 



19l|27|. 



'^Refs.rte, 
^Ref. 25 . 

'^Ref. [13. l—l 
'^The spin-dependent matrix elements from |l6l Il7 1 are used. 
•^Shell-model value with siv?9w = 0.23. 
■^This value was obtained by scaling the analytical result from 14t| km = 0.0049 by a factor 1.5. 
^Contains a 1.6% correction for finite nuclear size; the raw value is 0.094(16). 



Combining the experimental results for A [Im(£'pNc) / with our values for the 
spin-dependent matrix elements from Table EJ we obtain k = 0.117(16). The uncertainty 
comes from the uncertainty in the experimental value on the left side of Eq. fll4|l : the 
uncertainty in /? is negligibly small. In Ref. k = 0.112(16) was obtained by combi ning 
the same experimental data [l|, |26f with spin-dependent matrix elements from Refs. 
This value is also used in Refs. [isl. Iiol I27I. Differences with our value of k come only from 
differences in {7sFf\\z\\6sFj)^'^'°'\ 

Combining the effective value for k^ with the value K2 = 0.0140 from and the 
value K = 0.117(16) obtained above leads to Ka = 0.098(16), which is 8% larger than the 
value Ka = 0.090(16) obtained by Haxton and Wieman |lj] and 6% larger than the value 
Ka = 0.092(16) from Flambaum and Murray l28l. To clarify the sources of these differences, 
we compare our results with t 



15 



281 




hose from Refs. [1^ IM IJJ, |27|, |28|, |29| in Table |ITT1 We scaled 
to represent them in terms of the coefficients k, Ka, K2, 



the constants given in 
Khf used here. 

The revised value of k^ and {7sFf\\z\\6sFj)^'^'"''> obtained in this work increases the value 
of the ^^'^Cs anapole moment, and thereby slightly increases the differences between var- 
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ious experimental constraints on PNC meson coupling constants discussed in Since 
correlation corrections to {7sFf\\z\\6sFj)^'^'"'^ are large, 25% at RPA level, further accurate 
calculations are clearly desirable. 

Measuring the PNC electric-dipole transition between ground-state hyperfine levels is 
a potentially fruitful method for obtaining experimental anapole moments for atoms other 
than Cs. Schemes have been proposed to carry out such measurements and calculations have 



been carried out for various atoms in Refs. 



3G 



33, 



34l |. The contribution of the spin- 



independent interaction from Z exchange, which dominates the PNC dipole matrix element 
between different atomic levels, vanishes for the microwave transitions between hyperfine 
states of the same level. 

As an aid to the analysis of these microwave experiments, we give reduced matrix el- 
ements |z| induced by the spin-dependent interaction of Eq. (P) together with 
values of the third-order dipole matrix element |z| |-F/)'^'^^) for atoms of potential exper- 
imental interest in Table IIVI The corresponding calculations were carried out at the RPA 
level of approximation. The ratio of matrix elements again gives km- The ground state 
configurations of the atoms listed in the table are ns 1/2 or npi/2, and the hyperfine levels 
have angular momentum F = I ± 1/2. For some of the atoms considered in Table IIV|, 
RPA correlation corrections to the weak-hyperfine interference matrix element are no longer 
small; they contribute 20% and 34% to {Ff\\z\\Fj)''^^^ for Au and Ra+, respectively. 



IV. SUMMARY 



We have considered the PNC dipole matrix elements induced by the combined hyperfine- 
weak interaction and found that they are, at the few percent level of accuracy, proportional 
to the PNC dipole matrix elements induced by the spin-dependent interaction of Eq. (P), 
independent of Fj and Fp, for transitions Fj — Fp between hyperfine levels. The proportion- 
ality is not the result of an operator identity, but of similar angular momentum structures 
for the respective matrix elements. By carrying out calculations at the RPA level of ap- 
proximation, which are expected to be accurate to a few percent, we are able to extract 
effective coupling constants from the calculations. Although the dominant matrix ele- 
ment {wIFp\\z\\vIFf)^'^'°''' is sensitive to correlation corrections, increasing by 10-30% in Cs 
and Fr when correlation corrections are included, the matrix element {wIFf\\z\\vIFi)^^^'' is 
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TABLE IV: We list values of K^f for microwave transitions between ground-state hyperfine levels 
Fp — Fj in atoms of potential experimental interest. In this table, A is the atomic number, / is the 
nuclear spin, fij is the nuclear moment, Qw is the weak charge (including radiative corrections), 
C is the 50% fall off radius of both the nuclear p{r) and magnetization M(r) distributions (the 
10%-90% fall off distance is taken as 2.3 fm). The ground state configurations of the atoms 
considered here are nsx/2 or the hyperfine levels have angular momentum F = / it 1/2. 

The PNC reduced dipole matrix elements induced by the spin-dependent Hamiltonian of Eq. 
{Fp\\z\\Fj)^'^'^\ are given together with the third-order dipole matrix elements (Fpl l-^/)*-^^^; their 
ratio gives Khf. These calculations are carried out at the RPA level of approximation. Numbers in 
square brackets represent powers of 10. 



Atom 


A 


I 




Qw 


C(fm) 


nl 


Fi- 


Fp {Fp\\z\\Fi)M 


{Fp\\z\\Fi)i^') 




K 


39 


3/2 


0.39149 


-18.39 


3.611 


As 


1 - 


2 


-2.222[-13] 


-1.113[-16] 


5.01 


-04] 


K 


41 


3/2 


0.21448 


-20.36 


3.611 


As 


1 - 


2 


-2.222[-13] 


-6. 753 [-17] 


3.04 


-04] 


Rb 


85 


5/2 


1.3534 


-44.75 


4.871 


5s 


2 - 


3 


-2.550[-12] 


-5.432[-15] 


2.13 


-03] 


Rb 


87 


3/2 


2.7515 


-46.73 


4.871 


5s 


1 - 


2 


-1.363[-12] 


-1.027[-14] 


7.54 


-03] 


Cs 


133 


7/2 


2.5826 


-73.09 


5.675 


6s 


3 - 


4 


-1.724[-11] 


-7.791 [-14] 


4.52 


-03] 


Ba+ 


135 


3/2 


0.83863 


-74.01 


5.721 


6s 


1 - 


2 


-6.169[-12] 


-2.217[-14] 


3.59 


-03] 


Ba+ 


137 


3/2 


0.93735 


-75.98 


5.721 


6s 


1 - 


2 


-6.169[-12] 


-2.544[-14] 


4.12 


-03] 


Au 


197 


3/2 


0.14816 


-110.88 


6.554 


6s 


1 - 


2 


-1.601[-11] 


-1.912[-14] 


1.19 


-03] 


Tl 


203 


1/2 


1.6222 


-114.69 


6.618 


6pi/2 


- 


1 


-3.000[-ll] 


-3. 437 [-13] 


1.15 


-02] 


Tl 


205 


1/2 


1.6382 


-116.66 


6.618 


6Pl/2 


- 


1 


-3.000[-ll] 


-3.531[-13] 


1.18 


-02] 


Fr 


211 


9/2 


4.00 


-116.23 


6.733 


7s 


4 - 


5 


-2.379[-10] 


-2.223[-12] 


9.34 


-03] 


Fr 


223 


3/2 


1.17 


-128.08 


6.834 


7s 


1 - 


2 


-5.820[-ll] 


-5.187[-13] 


8.91 


-03] 


Ra+ 


223 


3/2 


0.2705 


-127.02 


6.866 


7s 


1 - 


2 


-5.987[-ll] 


-1.258[-13] 


2.10 


-03] 



correlation insensitive, changing by less than 6% for these cases. 

For the case of ^'^^Cs, the value of Khf is about 40% smaller than that obtained in an 



1 [l5| ar 



earlier calculation ll3 and slightly increases the size of the anapole moment of ^^^Cs inferred 
from experiment [l|, |l9|. Values of Khf are also presented for the 7s — 8s transition in Fr 
and for microwave transitions between ground-state hyperfine levels in atoms of potential 
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experimental interest. 
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APPENDIX A: THIRD-ORDER PERTURBATION THEORY 

We introduce a perturbation Vj = H^^^ + H^^^^ into the many-body Hamiltonian Hq 
describing an atom and expand the many-body wave function \1/ of the bound state v in 
powers of Vi 

^ = ^„ + + ^(2) ^ . . . 



to find 



*a) = 53MM^ (Ai) 



— En 



where 



\^ \^) MVil^) I ^ \m){m\Vi\n){n\Vi\v) 

^ {E,- E^f ^ 2^{E,- EJ{E, -EnY ^ ' 



E« = {v\Vi\v) = {v\H'^''''^\v) (A4) 



is the first-order correction to the energy. The approximate wave function + + ^>^v^ 
is normalized to second-order. 

The third-order matrix element of the dipole operator is given by 

{w\z\v)^''^ = {^^^M'^v) + (^iJ^kl^i^)) + )■ (A5) 
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Expanding this expression, we obtain 

^ {w\Vi\n) {n\Vi\m) {m\z\v) ^ ^--^ (w|V/|m)(m|z|n)(ri|V/|f ) 

"~ {Eu, — Em){Ew — En) ^ (E^ — En^iEy — En) 

v\Vi\n) {n\Vi\v) , , , \ i-i(i) ('"^I'^l'^)^'^!^^!'^) 



2^ {Ey-En)^ ^"^'"'"^^ 4- [Ey-En] 



2 



(-Et, — E„i) (Ey — 



^ \^ {w\z\n){n\Vj\m){m\Vi\v) 
+ 2^ _ R - ^^^^ 



Setting V/ = H^^^ + H^^^^ and retaining only those terms hnear in H^^\ we obtain the 
expression given in Eq. Q. It should be noted that the two terms above proportional to 
{w\z\v) do not contribute when the states v and w have the same parity. 

APPENDIX B: ANGULAR DECOMPOSITION 

The matrix element of the spin-independent operator H^^^ between single-particle states 
\i) and |j) is 

= ^^Qw 5.. dr [F,(r)G,(r) - G,(r)F,(r)] p(r). (Bl) 

Here, (k,, rrii) are angular momentum quantum numbers of the state \i) [/tj = =F(jj + 
1/2), for ji = ± 1/2], li and ji being the orbital and total angular momentum, respectively, 
of the state The functions Gi{r) and Fi{r) are the large and small radial components, 
respectively, of the Dirac wave function for the state We define the reduced matrix 
element of H^^'' as the coefficient of the angular momentum deltas in Eq. (jBlll . Using this 
(somewhat unconventional) definition, it follows that 

poo 

m^'^Wj) ='^Qwl dr mr)G,ir) - G,(r)F,(r)] p(r). (B2) 

We decompose the spin-dependent operators of the type H^''\ with k = (a, 2) in a 
spherical basis as^ 



^ In this equation, we omit the multiphcative factors Kq and K2 defined in the Sec.^to avoid confusion with 
the angular momentum quantum numbers Ki introduced in the previous paragraph. 
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The matrix element of the purely electronic operator Kj;!'^ between single-particle states \i) 
and |j) is 



G 



-{Kimi\a^j\ - K-jTUj) Gi{r)Fj{r) 



From this, it follows 



G 



{i\\K^'^)\\j)^i^ drp,{r) (-K,||a||«:,)F,(r)G,(r) 



-{Ki\\a\\-Kj) Gi{r)Fj{r)\^. 



Reduced matrix elements of the operator a are given by: 



1/2 1/2 li 

3j ji 1 
1/2 1/2 k 



(B3) 



(B4) 



(B5) 
(B6) 



In the above, we have used the notation I — and [j] — 2j + 1. 

The hyperfine operator is decomposed as 

X 

where /ix — Qi^nIx is the nuclear magnetic moment operator, and t\ is the electronic part 
of the hyperfine interaction. We may write the reduced matrix element of the magnetic 
moment operator in the nuclear ground state as 



{IMI) = Vn^ + l)(27 + l) g: pt^ . 



(B7) 



In the following, the factor //at is absorbed into the hyperfine interaction energy scale factor: 



lei |e|^ 1 K 

^ = 1.987131 X 10"^ a.u. 



47reo 2Mp cOg 

The electronic part of the hyperfine interaction for a point nucleus in these units is given by 



(B8) 
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where C^^^{^r) is a normalized vector spherical harmonic 35|, p. 210]. For a distributed 
nuclear magnetization M{r), Eq. ()B8|) becomes 



tx = -iV2 m(r), 



(B9) 



where the function m(r) is given by 



An 



m{r) = — I dss^ M{s) 
Jo 



dss^M{s) -i- / dss^M{s 



For the Fermi-type distribution given in Eq. (jH)), we find 



m(r, r < C) 



ar"^ f C — r\ a^r ( C — r 



a 

C-r 
a 



and 



m(r, r > C) = 1 — 
1 

where A/" is 



ar 



AT 



r-C 



1 + ^.^ + 6-^3 - 



In the previous three equations, 



(BIO) 



(Bll) 



\n-l 



The reduced matrix element of the hyperfine operator t is 



{j\\t\\i) = iK, + K,){-K,\\C4K,) / —{Gj{r)F,{r)+Fj{r)G,{r))m{r), (B12) 

Jo ^ 

where Ckq{'r) is a normalized spherical harmonic. Finally, we note that the reduced matrix 
element of the dipole operator z is 



\z\\t) = {k,\\C,\\k,) / drr{G,{r)G,{r) + F,{r)F,{r)) 
Jo 



(B13) 
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With the aid of the above expressions, the reduced third-order matrix element corre- 
sponding to Eq. (jH)) is found to be 



{wIF4z\\vIF,)^'''^ = (?V/(/ + l)(2/ + l)[F,][Fjx 

1 1 1 / / 1 

jj jw -^11 jj jv Fv 



' {w\\H^^\i) m\j) rntWv) ^ ^{w\\z\\i) {3\\t\\v) 



j+w 



F F 1 



/ / 1 



jv jj I \ I jj jw Fw 



'^ {w\\t\\j) m^'^v) , sr^ {w\\H^')\\i) mj) {3\\A\v) 



(B14) 



It is interesting to compare the interference term with the the second-order reduced matrix 
element of the dipole operator associated with the spin-dependent terms = 2,a, 

(m/F„||2||'t-/F.,>"'''' = V7(7+T)(27+T)[?y[?yx 



E 





' F^ F, r 








II 










Jj jv \ 



+ (-1) 



F F 1 



/ / 1 [ {w\\K^'^)\\j){j\\z\\v) 
V j I \ \ j w 



^w f j 



(B15) 



We find that the first term in this expression goes over to the first term in the interference 
term under the replacement 



{w\\z\\j){j\\K^'^^\v) ^ 

(u;||/7«||z) (zll^l 



9i 



v) , ^(•ii-lklli) {i||ff<"||j) (Jllilli.) 



e„ - e,; 



(ei- - ej) 



(B16) 
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and that a similar correspondence can be made for the second term. The completely different 
dependence on the intermediate state j on the two sides of the above expression explains 
the state-dependence of the coefficient K^f. 
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